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Introduction 

TIus  document  describes  the  installation  and  use  of 
the  ICCanvas  3D  visualization  program,  an  into*- 
active  graphics  tool  that  runs  in  the  \ficrosoft 
Windows  environment.  ICCanvas  was  developed  to 
simplify  the  study  and  communication  of  the  results 
of  complex  meteorological  models  and  other  3D 
simulations  by  exploiting  the  increasirig  gr^hics 
o^rabilities  of  nuCTocomputers.  The  program  allows 
you  to  use  your  personal  computer  to  view  level 
air&ces  on  structured  3D  dat^  much  of  vdiich  may 
have  been  produced  by  large  hydrocodes  running  on 
workstations  and  supercomputers.  The  ICCanvas 
con^ressed  binary  file  format  for  3D  data  is  up  to 
100  times  more  compact  than  a  sinq)le  text  repre¬ 
sentation  of  the  same  data,  meaning  that  form^ 
unwiddy  3D  data  sets  can  easily  be  stored  and 
transferred  using  inexpensive  removable  media  and 
ordiruuy  modems. 

ICCanvas  has  been  used  to  display  density  distribu¬ 
tions  of  dust  and  water  vapor  in  3D  simuiations  of 
nuclear  blast  clouds  and  to  reconstruct  the  surftce 
geometry  of  an  anthropomorphic  phantom’s  skde- 
ton  firom  computed  tomography  images.  Together 
with  sophisticated  3D  rendering  iqrplications, 
ICCam^  has  hdped  to  produce  Ugh-quality  ani¬ 
mated  views  of  the  development  of  a  blast  cloud 
featuring  footprint  projections  and  multiple  semi¬ 
transparent  surfaces.  These  animations  were  pro¬ 
duced  on  videotape  using  exactly  the  same  software 
and  techniques  currently  being  employed  to  create 
special  effects  for  several  commercial  tdevision 
programs. 
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Installation 


System  Requirements 

ICCanvas  requires  Microsoft  Windows  version  3. 1  or  later.  Recommended;  an  80386  (or 
better)  processor,  a  math  coprocessor,  at  least  4  m^abytes  of  memory,  a  hard  disk  with  at 
least  3  MB  of  five  space,  and  a  mouse. 

Hard  Disk  Procedure 

The  ICCanvas  program  package  consists  of  the  files  ICCANVAS.EXE,  ICUB.OLL,  and 
ICHELP.HLP.  To  install  ICCanvas,  copy  these  three  files  fivm  the  distribution  disk  to  any 
conveniem  place  on  your  hard  drive.  Hie  only  restriction  is  that  ICUB.OLL  and  ICHELP.HLP 
must  be  in  the  current  Windows  search  path  vdien  ICCanvas  is  launched.  The  easiest  way 
to  ensure  this  is  to  put  all  three  files  in  the  same  directory  on  your  hard  drive,  as  illustrated 
in  the  following  example. 


Assuming  your  hard  drive  is  C  and  the  ICCanvas  distribution  disk  is  in  drive  A,  execute 
these  commands  at  the  DOS  pronqit  to  copy  the  package  to  your  hard  drive: 


C:\>  ad  ieeanvM 
C:\>  od  ieeanvM 
C:\iccanvas>  copy  a:iceanvaa.«]e* 
C:\iccanvas>  copy  a:lelib.dll 
C:\iccanvas>  copy  a:ich«lp.hlp 


create  an  ICOmvas  directory 
make  it  the  current  directory 
copy  the  program 
copy  the  DLL 
copy  the  he^  file 


If  you  would  like  to  install  the  data  files  from  the  distribution  didc,  emer  these  DOS 
commands: 

C:\iccanvas>  md  data  create  a (kita directory 

C:\iccanvas>  xcopy  a:data\*.*  data  copy  the  data 

Of  course,  you  can  also  perform  the  installation  from  within  A^mdows  by  using  File 
Manager  to  create  directories  and  copy  files.  To  make  the  ICCanvas  icon  visible  and 
selectable  in  Program  Manager,  create  a  new  program  item  (or  group)  by  selecting  New  in 
the  File  menu  of  the  Program  Manager  window.  If  you  copied  the  data  from  the 
distribution  disk,  you  will  probably  want  to  make  the  data  directory  the  ICCanvas  working 
directory. 


Defining  TEMP 

ICCanvas  creates  temporary  files  during  the  construction  of  a  surftce.  These  files  are 
used  to  collect  the  vertices  and  faces  that  make  up  a  surfiwe  while  the  data  file  is  being 
scanned.  Once  the  surface  is  built,  these  files  are  no  longer  needed,  and  ICCanvas  deletes 
them.  ICCanvas  asks  Windows  which  disk  should  be  used  for  temporary  files.  Windows 
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will  oidiiurily  choose  the  first  hard  drive  it  finds  on  your  system,  but  you  may  tell 
Windows  to  use  another  drive  by  defining  a  TEMP  environment  variable.  (Note, 
however,  that  ICCanvas  does  not  require  ymi  to  do  this.  The  procedure  described  in  the 
foDowing  paragraphs  may  help  you  improve  the  performance  of  ICCanvas.  If  you  aren’t 
comfortable  with  DOS  or  with  the  techniques  described,  you  may  safely  ignore  this  part.) 

Ymi  set  an  environment  variable  by  using  the  DOS  coimnand  Mt.  This  command  can  be 
entered  at  the  DOS  pron:q>t,  but  it  usually  appears  as  a  line  in  your  autCMEXEC.bat  file,  the 
batch  file  that  is  executed  every  time  you  turn  on  your  computer.  If  you  have  set  up  a 
RAM  drive,  you  can  improve  the  performance  of  ICCanvas  by  making  the  TEMP  variable 
point  to  your  RAM  drive.  Assuming  the  RAM  drive  is  D,  the  Mt  command  would  be 

Mt  TEMP^:\ 

If  you  insert  this  command  in  autc£XEC.Bat,  Wmiows  will  tell  ICCanvas  to  use  your 
RAM  drive  for  hs  temporary  files. 

There  are  a  number  of  issues  you  should  be  aware  of  before  setting  a  TEMP  environment 
variable.  Many  Windows  programs  that  use  temporary  files  either  look  for  TEMP 
themselves  or  ask  Windows  to  do  h,  and  in  Standard  mode,  Wndows  itself  may  use  the 
TEMP  drive  for  iqiplication  swap  files  if  the  WIN.INI  variable  •wapfll*  is  undefined. 
Your  TEMP  drive  must  be  large  enough  to  accommodate  the  temporary  files  fi’om  each  of 
these  programs,  and  this  can  be  of  some  concern  if  the  drive  is  a  Sxed-aze  RAM  drive, 
especially  since  larger  RAM  drives  reduce  the  amoum  of  firee  memory  available  to 
Windows. 

In  practice,  we’ve  found  that  a  512K  RAM  drive  woilcs  wdl  and  is  entL%  adequate  for 
most  applications.  ICCanvas  requires  a  maximum  of 270K  for  its  temporary  files,  and  this 
space  is  only  required  while  ICCanvas  is  actually  building  a  surfiu:e.  Interactions  between 
the  TEMP  environment  variable,  RAM  drives,  swap  files,  drive  caching  and  Windows 
memory  management  are  discussed  at  some  length  in  your  Windows  documentation 
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Execution  Procedures 

ICOuivas  is  an  event’^bfven  program.  It  acts  only  in  response  to  requests  from  the  user. 
As  those  with  experience  tunning  Windows  programs  know,  this  makes  ICCanvas  more 
flexible  than  an  entire^  sequential  program  but  puts  greater  responsibility  on  the  user,  on 
whrnn  fiJls  the  task  of  tdl^  the  program  wduh  to  do.  This  section  therefore  begins  by 
trying  to  give  you  an  idea  of  what  the  program  can  do.  It  also  lists  and  describes  each 
option  ICCanvas  offers  you  >^e  the  program  is  running. 

Betting  Started 

ICCanvas  allows  you  to  visualize  three-dinwnsional  data.  It  assumes  that  the  data 
represents  a  3D  array  of  density  values  that  model  a  volume.  After  examining  tins  data,  it 
is  able  to  build  surfrices  that  form  boundaries  around  higher  density  parts  of  the  volume. 
You  control  the  density  value  at  ^ch  the  surfrice  is  buih,  along  with  your  point  of  view 
and  the  method  us^  to  draw  the  surface. 

To  begin,  use  the  Open 
menu  option  in  the  File 
menu  to  open  a  data  file. 

ICCanvas  examines  the 
file  you  choose  and  shows 
you  an  approximate  dis¬ 
tribution  of  the  densities 
in  the  file  using  the  Level 
dialog  By  default,  the 
level  is  set  to  the  median 
value  in  this  distribution, 
but  you  may  set  it  to  what 
you’d  like.  Densities 
higher  than  the  level  you 
choose  will  be  inside  the 
resulting  surface. 

Once  the  data  file  is  chosen  and  the  level  set,  ICCanvas  reads  the  data  and  constructs  a 
surface.  The  progress  of  the  build  is  displayed,  and  during  most  of  that  time  you  have  the 
option  of  canceling  the  build.  Assuming  you  don’t  cancel,  the  surftce  is  then  drawn. 

To  view  the  surface  from  different  angles,  click  and  hold  the  left  mouse  button  anywhere 
in  the  ICCanvas  display.  A  picture  of  a  square  with  a  vertical  line  connected  at  one  comer 
(the  axis  image)  will  be  drawn  on  top  of  the  sur&ce.  Moving  the  mouse  will  change  the 
orientation  of  the  axis  image,  and  when  you  release  the  left  mouse  button  the  surfrice  will 
be  redrawn  from  the  new  viewing  an^e.  For  greater  control  over  the  view,  use  the 
Camera  and  Center  menu  options  in  the  Settings  menu. 
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To  build  a  suiftce  fiom  the  same  data  but  with  a  different  level,  use  the  Level  menu 
option  in  the  Settings  menu.  To  buUd  a  surface  from  new  data,  just  Open  the  new  data 
file.  Note  that,  unless  you  change  it,  the  level  setting  will  remain  the  same  as  it  was  for  the 
dd  data  file,  and  the  old  level  may  not  be  ^ipropriate  for  the  new  data. 

The  surfiice  can  be  drawn,  or  rendered,  in  one  of  three  ways  controlled  by  the  Render 
Style  menu  options  in  the  Settings  menu.  The  ICCanvas  display  can  be  saved  to  a  file  in 
one  of  three  images  formats  using  the  Save  Image  menu  options  in  the  File  menu.  The 
display  can  also  be  posted  to  the  Windows  clif^ard  with  the  Copy  Image  menu  option  in 
the  E^t  menu.  The  actual  geometry  of  the  sur&ce  can  be  saved  in  either  of  two  formats 
using  the  Save  Surface  menu  options  in  the  File  menu.  Using  the  Window  Size  and 
Colon  menu  options  in  the  Settings  menu,  you  can  predsely  set  the  size  of  the  display  and 
modify  the  colors  of  the  background  and  the  ‘‘floor,**  a  rectangle  that  provides  a  visual  cue 
to  the  orientation  of  the  surface. 

The  Help  menu  gives  you  access  to  much  of  the  text  of  this  section  online,  as  well  as  an 
About  display  containing  a  copyright  notice  and  acknowledgments  of  support  for  the 
development  of  the  program. 

Once  you’ve  finished  with  ICCanvas,  choose  Exit  fiom  the  File  menu  to  quh. 


Open 

You  use  the  Open  menu  option  to  tell  the  program  the  name  of  a  new  IC  data  file.  This  is 
the  first  step  in  constructing  a  surface  from  new  data.  For  historical  reasons,  IC  format 
data  files  on  the  distribution  disk  all  end  with  extensions  of  .  RLE,  but  if  you  happen  to 
have  IC  format  files  with  different  extensions,  you  can  see  them  in  the  Open  dialog  by 
asking  for  All  in  the  List  Files  of  Type  drop>down  list.  After  you  specify  a  data  file  and  set 
the  density  level  in  the  Level  dialog,  the  data  file  will  be  read  by  the  program  to  find  and 
build  a  surface  of  constant  density  through  the  data. 

If  you  choose  the  Cancel  button  in  the  Open  dialog,  the  data  file  you  were  working  with 
previously  will  remain  the  current  file.  If  you  choose  OK,  the  current  data  file  is  forgotten 
and  its  place  taken  by  the  new  one  you  specify.  Before  the  program  builds  a  sur&ce, 
you’ll  be  given  a  chance  to  change  the  current  level.  The  levels  listed  in  the  Level  dialog 
will  reflect  the  distribution  in  the  new  file. 

Each  IC  data  file  contains  a  3D  array  of  den«ty  values  in  a  compressed  binary  format, 
along  with  information  about  the  distance  between  data  points,  the  origin  of  the  data’s 
coordinate  system,  the  distribution  of  densities,  the  date  and  time  when  the  data  was 
created,  and  if  appropriate,  the  time  point  the  data  represents.  Si  s  Appendix  B,  IC  File 
Format,  for  more  information. 
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Save  Image 

You  use  the  Save  Image  menu  options  to  create  a  file  containing  the  image  currently  in  the 
program’s  window.  The  image  can  that  be  loaded  into  other  programs  for  editing, 
presentation  or  printing,  or  can  be  given  to  someone  else.  The  size  of  the  image  and  the 
number  of  colors  it  contains  are  determined  by  your  display  hardware  and  the  Windows 
video  driver  you  are  using.  Images  can  be  saved  in  one  of  three  formats,  each  of  which 
has  a  broad  base  of  support  among  gn^rhics  replications  for  personal  computers. 

Windows  BMP.  Files  of  this  type  contain  Device-Independent  Bitmaps,  a  format 
introduced  in  Windows  3.0  and  based  on  the  bitmap  format  used  in  OS/2  Presentation 
Manager  1.1.  Some  support  for  BMPs  is  buih  into  Windows,  so  many  Windows 
programs  can  read  and  write  them.  ICCanvas  writes  uncompressed  BMPs. 

PC  Paintbrush  PCX.  This  is  a  format  creat  xi  for  ZSofi’s  PC  Paintbrush  and  widely 
supported  by  both  DOS  and  Windows  programs.  ICCanvas  writes  PCX  version  S,  which, 
unlike  earlier  versions,  aUows  8-  and  24-bit  images.  Because  PCX  uses  compression,  this 
format  produces  files  that  are  almost  always  sgnificantly  smaller  than  BMPs. 

Macintosh  PICT.  The  PICT  format  is  built  imo  the  Apple  Macintosh  operating  system. 
ICCanvas  writes  pixels  in  0x(X)9A  (24-bit)  and  0x0098  opcode  blocks,  both  of  which  are 
compressed. 


Save  Surface 

You  use  the  Save  Surface  menu  options  to  create  a  file  containing  the  geometry  of  the 
current  surface.  ICCanvas  defines  a  surface  as  a  mesh  of  triangles,  each  of  which  is 
determined  by  the  (x,  y,  z)  positions  of  its  three  vertices.  Unlike  the  Save  Image  options. 
Save  Surface  preserves  all  of  the  3D  information  needed  to  reconstruct  the  surface  and 
view  it  from  any  angle.  These  geometry  files  can  be  loaded  into  CAD  programs  or 
programs  that  perform  photorealistic  rendoing,  scioitific  visualization,  or  finite  element 
analysis. 

AutoCAD  DXF.  This  is  the  format  used  by  3D  applications  from  Autodesk  and  widely 
supported  by  other  CAD  and  3D  rendering  programs.  Its  wide  support,  in  &ct,  is  the 
primary  reason  for  its  inclusion  in  ICCanvas.  You  may  find  that  in  other  respects 
(unwieldy  file  sizes  for  typical  IC  surfiices,  for  instance)  DXF  isn’t  entirely  satisfiictory. 

Wavefront  OBJ.  OBJ  is  the  text  version  of  geometry  files  used  by  Wavefiont 
Technologies’  Advanced  Visualizer  and  otho'  programs.  Support  for  this  format  will  be 
found  predominantly  on  Silicon  Grqihics  workstations  and  in  some  3D  modeling  and 
rendering  programs  for  the  Apple  Macimosh,  but  programmers  vidio  write  thdr  own 
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^)plicatioiis  to  read  ICCanvas  geometry  files  will  prefer  the  simplicity  and  relative  speed 
of  this  format  to  DXF. 

Coiqf  Image  lilllllillLl 

You  use  the  Copy  Image  memi  option  to  put  a  copy  of  the  image  currently  in  the  display 
into  the  Windows  dipboard.  You  can  then  switch  to  a  diSerem  program  and  use  that 
program’s  Paste  fimcticn.  This  is  an  easy  way  to  put  ICCanvas  graphics  imo  a  paint 
program,  a  presentation,  or  a  word  processing  document. 

Note  that  this  will  woric  even  if  you  quit  ICCanvas  before  performing  the  paste.  ICCanvas 
puts  a  permanent  copy  of  the  image  into  the  clipboard;  it  has  a  life  of  its  own  and  doesn’t 
require  ICCanvas  for  rendering. 

As  is  customary  with  the  clipboard,  ICCanvas  ddetes  the  previous  clipboard  contents 
before  posting  hs  itiukge.  This  means  you  should  only  Copy  if  you  no  longer  need  what  is 
already  there,  including  previous  ICCanvas  images. 


Camera 


]gfla4awSlM_ 
Cgliw _ » 


You  use  the  Camera  dialog  to  control  the  position  and  orientation  of  the  ICCanvas 
camera.  The  ’’camera”  is  just  a  convoiient  metaphor  that  disguises  many  of  the  details  of 
the  viewing  system’s  geometry,  a  collection  of  vectors  and  transformations  that 
determines  what  gets  drawn  in  the  display. 

Azimuth,  Latitude.  These  are  used  to  move  the 
camera  around  the  current  surface.  A  change  in 
azimuth  is  equivalent  to  east-west  motion  around  a 
globe,  while  a  change  in  latitude  is  equival«it  to 
north-south  motion  The  axis  im^e  (the  picture 
of  a  square  patch  of  ground  with  a  vertical  a»s 
attached  at  the  origin)  will  change  to  reflect  azimuth  and  latitude  values  you  emer. 
Conversely,  moving  the  axis  image  with  the  mouse  will  change  the  values  in  the  text 
boxes 


Zoom.  As  the  name  implies,  this  setting  allows  you  to  magnify  or  reduce  the  image  in  the 
display  Unlike  a  real  camera,  this  setting  has  no  effect  on  the  amoum  of  perspective  in 
the  image. 

Prcseu.  Six  buttons  are  provided  for  special  camera  settings  that  align  the  view  with  an 
axis,  set  the  Zoom  to  1 .0,  and  set  the  Center  to  the  center  of  the  data  space.  These  are 
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useful  as  standard  points  of  view,  and  may  also  be  hdpfiil  after  a  long  series  of  camera 
maniinilations  have  left  you  a  little  uncertain  about  whm  you  are. 

The  camera  dialog  is  available  r^ardless  of  whether  you  have  a  sur&ce  to  be  displayed.  If 
tlwre  is  a  current  sur&ce,  the  text  of  the  OK  button  will  change  to  Render,  a  reminder  that 
accq>ting  the  changes  you  make  will  cause  the  display  to  be  redrawn  using  the  new 
settings.  While  the  Camera  dialog  is  displayed,  you  may  press  the  Reset  button  at  any 
time  to  restore  the  settings  to  vduit  they  were  v^en  the  dialog  was  first  opened. 


The  azimuth  and  latitude  settings  can  be  changed  without  using  the  Camera  dialog.  If  you 
click  and  hold  the  left  mouse  button  anywhere  in  the  ICCanvas  di^lay,  the  axis  image  will 
tq)pear.  For  as  long  as  the  left  mouse  button  is  held  down,  you  may  move  the  axis  image 
by  moving  the  mouse.  You  accept  the  new  camera  orientation  when  you  release  the  left 
mouse  button. 


Level 


CjlOT _ » 


You  use  the  Level  dialog  to  specify  the  density  at  wduch  to  build  an  isosur&ce.  Once  a 
level  is  defined,  ICCanvas  will  scan  the  currem  data  file  (the  one  you  most  recently 
opened)  to  find  the  path  of  a  surface  of  constam  density  at  that  level.  The  volume 
enclosed  by  this  surface  contains  all  the  data  points  with  values  at  or  above  the  level  you 
choose. 

ICCanvas  provides  a  table  of  levels  that  tries  to  pve  some  indi* 
cation  of  the  distribution  of  denaties  in  the  data.  Each  level  is 
paired  with  a  percentage  that  tells  how  many  of  the  non-zero 
data  points  will  be  enclosed  by  a  surfi^e  at  that  level.  In  this 
respect,  ICCanvas  assumes  that  the  data  represems  a  mass  of 
some  arbitrary  shape  sitting  in  a  rectangular  void.  The  statis¬ 
tics  in  the  list  attempt  to  be  relevant  to  the  mass  and  not  to  the 
void. 


You  can  select  one  of  the  listed  levels  or  type  one  of  your  own.  The  value  you  emer  must 
be  an  actual  density  value  and  not  a  percentage.  It  may  be  entered  either  in  the  power-of- 
10  notation  used  in  the  list  or  as  an  ordinary  decimal  number,  but  for  most  purposes  it 
should  be  in  the  range  offered  by  the  list,  since  this  indicates  what  values  are  in  the  data. 
Note  that  levels  are  expressed  in  arbitrary  units  that  depend  on  the  data  file. 

The  level  you  set  will  remain  in  effect  until  the  next  time  you  set  a  level.  You  may  find 
that  a  level  you  enter  will  be  slightly  altered  the  next  time  you  see  h.  ICCanvas  works 
with  densities  as  5-place  common  logarithms,  and  this  will  result  in  round-oflf  errors  that 
are  harmless  but  may  occasionally  look  peculiar.  For  a  more  technical  explanation,  see  the 
IC  File  Format  topic. 
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You  use  the  Center  menu  option  to  center  the  display  around  a  point  you  choose.  When 
you  activate  centering,  the  mouse  cursor  changes  to  a  plus  image  a^  a  small  window 
opens  to  show  you  where  you  are  in  object  coordinates.  To  move  the  cemer  of  the 
display,  move  tlw  cursor  to  the  new  center  poim  and  click  the  left  mouse  button.  To 
cancel  centerii^  with  no  change,  click  the  right  nK)use  button. 

Imemally,  ICCanvas  treats  the  cemer  as  a  3D  poim  used  to  position  the  camera.  The 
center  is  the  origin  of  the  azimuth  and  latitude  rotations,  or  the  pivot  point  around  which 
the  camera  moves.  You  can  get  complete  control  of  the  3D  position  of  the  cemer  poim  by 
u^g  Center  in  combination  with  the  preset  orientations  in  the  Camera  dialog,  in  the 
following  way; 


1 .  In  the  Camera  dialog,  choose  the  Top  or  Under  view  to  get  an  x-y  projection. 

2.  From  the  Settings  menu,  choose  Cemer  to  set  the  x  and  y  coordinates. 

3.  In  the  Camera  dialog,  change  the  latitude  to  0**  to  get  a  projection  that  includes  z. 
Leave  the  atimuth  at  O'*. 

4.  Choose  Center  agun  to  set  the  z  coordinate. 

You  should  avoid  using  the  preset  orientations  once  you  have  a  cemer  that  you  want  to 
keep.  The  presets  automatically  restore  the  cemer  poim  to  the  center  of  object  space, 
which  is  why  they  weren’t  used  in  stq}  (3)  id>ove. 


BEB _ 

1 

1 

SinWii. 

Ccfliir-. 

1 

Window  Size 

1  Cfian  >1 

The  Window  Size  dialog  allows  you  to  predsely  set  the  size  of  the  ICCanvas  display. 
This  is  useful  for  controlling  the  size  of  the  images  created  by  the  Save  Image  and  Copy 
menu  options.  Window  Size  controls  the  dimensions  of  the  display  area  (in  Windows 
parlance,  the  client  area  of  the  window),  and  not  the  entire  window. 

You  can  set  the  size  in  either  pixels  (picture  dots)  or  inches. 

Sizes  expressed  in  inches  depend  on  a  scale  factor,  the  number  of 
pixels  per  inch.  The  video  driver  for  your  di^lay  provides  a 
default  scale  that  is  displayed  in  the  DPI  (dots  per  inch)  text  box 
the  first  time  the  Wmdow  Size  dialog  is  open^.  Changing  the 
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DPI  setting  afiects  the  scaling  ICCanvas  performs  to  determine  pixd  dimensions  from 
sizes  e^qiressed  in  inches. 

Suppose  you  would  Uke  to  produce  a  3"  x  3*  ICCanvas  image  for  inclusion  in  a  word 
processing  document.  Your  word  processor  allows  you  to  set  the  DPI  for  an  image,  and 
you  know  you  can  get  good  results  from  your  printer  if  images  are  printed  at  ISO  DPI.  In 
ICCanvas,  you  can  set  the  DPI  to  150  and  the  window  dimensions  to  3"  x  3",  and  the 
display  wffl  be  sized  to  the  correct  pixel  dimensions. 

Note  that  the  DPI  setting  is  internal  to  ICCanvas.  It  has  no  effect  on  Windows,  nor  is  it 
explicitly  written  imo  image  files  you  save  at  images  you  copy  to  the  dipboard. 


Render  Style  ‘ 


You  use  the  Render  Style  menu  option  to  set  the  way  the  surfoce  is  drawn  in  the 
ICCanvas  display.  The  chdoe  involves  a  tradeoff  between  the  quality  of  the  image,  the 
information  it  contains,  and  the  time  h  takes  to  draw.  The  new  render  style  will  be 
employed  the  next  time  ICCanvas  draws  the  surfiwe.  To  see  the  change  immediately,  click 
the  left  mouse  button  anywdiere  in  the  ICCanvas  diq>lay. 


Facet  This  option  uses  constant  shadii^  aat>ss  each  of  the  triangles  that  make  up  the 
surface.  The  rendering  is  relativdy  fiut,  but  the  sinq>le  shading  enqrhasizes  the  polygonal 
nature  of  the  surfoce  iq>proximation,  which  may  dishact  you  fit>m  features  of  the  siufece 
that  should  ^pear  continuous.  On  displays 
with  16  or  fewer  colors,  fecet  rendering  uses 
the  ordered  dither  method  of  color  approx' 
imation  built  into  Windows,  and  this  may 
also  be  distracting. 

Smooth.  This  option  employs  Gouraud 
interpolated  shading  to  draw  the  sur&:e 
triangles.  Interpolated  shading  ensures  that 
the  shading  of  the  surface  is  continuous 
across  triangle  edges.  If  the  polygon  mesh  is 
an  approximation  of  a  smooth  surfece,  con> 
tinuous  shading  across  polygon  boundaries 
provides  a  better  representation  of  the  un¬ 
derlying  surface.  Gouraud  shading  works  at  the  level  of  jncture  dots,  or  pixels,  rather  than 
triangles,  so  it  can  be  somewhat  slower.  In  ICCanvas,  h  also  uses  absolute  color,  and  on 
displays  with  16  or  fewer  colors  this  can  mean  that  the  surfece  is  drawn  with  only  the 
three  or  four  shades  of  gray  that  happen  to  be  available. 
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Wire.  This  opdon  draws  a  wirefiame  version  of  the  surftce.  Wirefiames  wVf  the 
polyg<Mi  representation  of  the  surftce  explicit  and  may  be  hdpfiil  in  clarifying  the  gf 
certain  feahires  in  cases  t^iiere  the  shaded  renderings  ^>pear  amhig^ioys  xhis  is  also  the 
fittest  renderii^  method. 

The  shadii^  ofthe  surftce  for  the  Facet  and  Smooth  options  uses  the  Lambert  relation  for 
difiiise  reflected  light,  in  which  the  apparent  brightm^  of  a  fiu:e  is  proportional  to  the 
cosine  of  tlw  angle  formed  by  the  sur&ce  normal  and  vector  pointing  toward  the  light 
source.  In  ICCanvas,  the  light  source  is  always  coincident  with  the  camera. 


Colors 

You  use  the  Colors  menu  options  to  set  the  colors  of  the  background  and  the  floor  in  the 
ICCanvas  display.  This  is  useful  if  you  plan  to  print  ICCanvas  images  in  black  and  white 
and  would  like  a  white  background,  for  instance.  You  may  also  find  that  the  defiuih 
ICCanvas  colors  are  not  good  choices  for  your  particular  di^lay,  or  that  th^  are  just 
unappealing  to  you. 

Color  changes  remain  in  effect  only  for  the  current  ICCanvas  session. 
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SECTION  4 

Error  Messages 

This  section  lists  the  enx>r  messages  you  might  see  while  using  ICCanvas.  The  messages 
are  listed  in  execution  order — messages  for  errors  that  would  occur  early  in  an  ICCanvas 
session  are  listed  first. 

•  File  Erron  Cannot  find  ICLIB.DLL 

•  File  ErroR  Cannot  find  COMMDLG.DLL 

When  Wmdows  starts  ICCanvas,  it  also  tries  to  open  ICLIB.DLL  and  COMMDLG.DLL, 
since  ICCanvas  requires  both  files  in  order  to  run.  ICLIB.DLL  is  the  1C  dynamic-link 
library  and  (X>ntains  much  of  the  fiinctionality  of  the  program.  COMMDLG.DLL  is  pan 
of  the  Windows  3.1  package  and  contains  the  Open  and  Choose  Color  common  dialogs 
Both  files  must  be  in  the  Wmdows  or  WmdowsXSystem  directories,  or  in  the  same 
directory  as  ICCANVAS.EXE,  when  you  suut  ICCanvas.  See  the  Installation  section  for 
more  information.  ‘ 

•  Filename 

Cannot  find  this  file. 

Please  verify  that  the  correct  path  and  filename  are  given. 

The  filename  you  entered  in  the  Open  dialog  can’t  be  matched  to  an  existing  file.  If  the 
problem  with  the  name  of  the  data  file  isn’t  obvious,  try  entering  the  name  by  picking  the 
path  and  file  from  the  lists  provided  in  the  Open  dialog,  rather  than  typing  it  into  the 
filename  edit  field. 

•  Unable  to  open  "filename*^. 

Windows  thinks  that  the  data  file  you  chose  fiom  the  Open  dialog  exists,  but  ICCanvas  is 
unable  to  open  it.  This  won’t  happen  often,  although  h’s  possible  that  read  access  to  a  file 
might  be  denied  because,  for  example,  another  application  has  the  file  open. 

•  Filename 

is  not  an  1C  format  file. 

The  file  you  chose  from  the  Open  dialog  is  not  a  l^al  data  file.  ICCanvas  checks  the  file 
you  name  to  ensure  that  it  is  large  enough  to  be  a  data  file,  that  the  name  field  is  a  legal  C- 
language  string,  and  that  the  dimensions  are  reasonable.  This  will  prevent  you,  in  most 
cases,  from  accidentally  loading  a  non-IC  file  as  a  data  file. 

•  General  Protection  Fault 

Windows  displays  this  message  when  an  application  has  attempted  to  read  firom  or  write 
to  memory  it  doesn’t  own.  In  ICCanvas,  this  is  usually  a  symptom  of  a  bad  data  file  (a 
damaged  or  non-IC  file),  especially  if  it  occurs  during  the  construction  of  a  sur&ce.  In 
any  case,  this  is  a  serious  message,  and  Avhhout  knowing  the  cause,  you  shouldn’t  assume 
that  any  part  of  Windows  will  fiinction  reliably  once  you  have  received  it — ^while  protected 
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mode  prevents  most  kinds  of  illegal  addresang,  this  is  only  protection  against  a  symptom 
and  not  whatever  pathology  caused  it.  It  is  safest  to  exit  and  restart  Windows. 

•  The  surface  is  too  large  at  the  specified  Level 
Try  using  a  higher  (less  inclusive)  value. 

Surfaces  cannot  contain  more  than  32,767  triangles  and  16,384  points,  a  practical  limit 
that  arises  both  because  of  the  way  PCs  use  memory  and  as  a  judgmern  about  what 
ICCanvas  can  reasonabfy  do.  This  message  may  also  result  from  the  use  of  an 
inappropriate  level  left  over  from  a  previous  sur&ce.  Guided  by  the  denaty  distribution 
di^layed  in  the  Level  dialog,  change  the  level  to  one  that  will  create  a  sur&ce  with  fewer 
cninp«n^s  and  try  the  build  ^ain.  Note  that,  until  the  build  is  successful  ICCanvas  has 
no  sur&ce  to  display. 

•  No  surface  exists  at  the  specified  Level 
Try  using  a  lower  (more  inclusive)  value. 

The  build  produced  no  triangles,  in^lying  that  all  of  the  data  values  are  less  than  the 
current  level.  This  message  may  re^  from  the  use  of  an  inappropriate  level  left  over 
from  a  previous  surface,  (juided  by  the  den^  distribution  displayed  in  the  Level  dialog, 
change  the  level  to  one  that  will  firid  a  surftce  and  try  the  build  again.  Note  that,  until  a 
non-empty  sur^e  is  built,  ICCanvas  has  no  surftce  to  display. 

•  ICLIB:  icBuildInitO  Opening  or  Reading  source. 

•  ICLIB:  icBuildInitO  Opening  temp  files. 

•  ICLIB:  icBuildInitO  Memory  allocation. 

•  ICUB:  icBuildInitO  Writing  vertices. 

•  ICLIB:  icBuildO  Reading  source  or  Writing  vertices  or  faces. 

•  ICLIB:  icCompactFaccsO  Memory  attocation  or  file  i/o. 

•  ICLIB:  icBuildFinishO  Memory  allocation. 

•  ICLIB:  icBuildFinishO  Reading  vertices  or  faces. 

These  errors  can  occur  during  the  construction  of  a  surface,  and  all  indicate  a  problem 
with  a  disk  or  memory  resource.  A  memory  allocation  error  usually  means  that  you  have 
run  out  of  free  memory,  but  it  may  also  be  caused  by  memory  fragmentati<m — ^Windows 
tries  very  hard  to  rearrange  items  in  memory  so  that  the  remaining  five  memory  is  in  a 
small  number  of  large,  contiguous  blocks,  and  for  various  reasons  h  may  not  always 
succeed.  ICCanvas  can  require  as  much  as  l.S  MB  (megabytes)  of  tnonory,  in  separate 
pieces  that  are  never  larger  than  64K  (kilobytes)  each.  If  you  are  running  other  programs 
concurrently  with  ICCanvas,  try  quitting  from  those  programs  and  restarting  the  build.  As 
a  last  resort,  you  can  check  for  memory  fragmentation  by  quitting  and  restarting  Windows 
itself,  but  be  aware  that  a  program  in  your  Windows  startup— a  screen  saver,  for 
example — might  be  the  source  of  the  problem. 

Disk  errors  are  less  common  and  can  be  caused  by  a  loss  of  access  to  a  disk  or  file,  a 
iiill  disk,  or  by  a  hardware  problem.  Loss  of  access  can  occur  when  files  are  being 
from  or  written  to  remov^Ie  media  (floppy  disks,  hard  disk  cartridges  or  tape,  for 
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insUnce)  or  storage  accessed  over  a  network.  ICCanvas  may  also  be  attempting  to  create 
hs  vertex  and  fiice  temporary  files  on  a  write-protected  or  nearly  full  disk.  During  the 
build,  ICCanvas  uses  up  to  270K  of  disk  space  for  its  tenq)oraiy  files.  These  files  are 
ddeted  when  the  build  is  complete. 

The  disk  used  for  temporary  files  is  determined  by  Windows,  but  you  may  tell  Windows  to 
use  a  particular  drive  fi^r  temporary  files  by  defining  a  TEMP  environmem  variable.  See 
the  Inkallation  section.  Note  to  programmers;  the  IC  DLL  rdies  on  the  API  function 
GetTempFileNameO  and  does  not  attempt  to  validate  the  name  or  check  disk  free  space. 

•  Unable  to  save  the  image. 

•  Unable  to  save  the  surface. 

These  messages  can  i4)pear  during  processing  of  the  Save  Image  and  Save  Surface  menu 
options.  The  errors  tto  cause  them  are  of  the  same  type  as  those  that  may  occur  during 
surface  construction:  a  lack  of  memory  or  a  problem  with  the  output  file.  In  particular, 
make  sure  that  the  Rename  you  supply  for  the  save  operation  is  valid — the  disk  exists,  it's 
inserted,  and  h  isn’t  write-protected,  the  path  exists  and  is  typed  correctly,  the  base  name 
isn’t  more  than  8  characters  long,  the  extenrion  isn’t  more  than  3  characters  long,  all  the 
characters  are  legal  for  filenames,  ymi  aren’t  trying  to  overwrite  an  existing  file  with  its 
write/delete  protection  bit  set.  etc.  yTindows  slxHild  alert  you  to  some  of  these  problems 
before  the  program  has  gotten  far  enough  to  issue  ttese  errors,  but  make  no  assungrtions. 

Also  make  sure  that  the  disk  has  enough  fiee  space  for  the  file.  Surfime  files,  in  particular, 
can  be  quite  large,  requiring  more  space  than  is  available  on  an  enqrty  high-density  floppy 
in  some  cases. 
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AlfBfflllA 

The  1C  Approach 

This  ^)pendix  gives  a  brief  devdopment  of  level  suiftces  and  discusses  the  method  IC 
uses  to  build  ismuiftces. 

ICCanvas  and  the  IC  DLL  construct  vibat  are  called  level  surfixes  or  isosurfaces.  For 
any  number  c,  the  set  of  points  {x,y,x)  for  which  /{x,y,x)  =  c  is  called  a  level  surface  of 
f .  If  f{x^y^x)  gives  the  temperature  at  pmnts  in  space,  the  level  surface 

f{x,y,x)^c  b  the  set  of  ail  points  at  vidnch  the  tenqierature  b  c  and  b  called  an 
isothermal  suiftce.  Similarly,  y(x,y,x)  ^  c  mi^  be  a  surftce  of  constam  voltage,  called 
an  equipotential  sui&ce. 

Spheres  centered  at  the  origin  are  a  set  of  levd  sui&ces  ofthe  fiinction 
/{x,y,x)  =  x^+y^+x* 

since  +y^  =  c  b  an  equation  of  sudi  a  sphere  for  non-zero  c.  In  thb  example,  the 

isosurface  can  be  constructed  analytically  becmse  the  fiinction  b  knowiL  IC,  on  the  other 
hand,  works  with  function  values  tabulated  at  r^ular  intervab  and  arranged  in  a 

3D  array,  or  scalar  field.  See  the  IC  File  Format  ^ipendix  for  infiirmation  on  the  way 
scalar  fields  are  presented  to  IC  as  data  files. 

IC  extracts  surftces  by  scanning  a  scalar  fidd  two  slices  at  a  time,  firom  the  bottom  up. 
Each  adjacent  pair  of  slices  forms  a  lattice,  with  data  values  where  the  grid  lines  meet  (or, 
the  other  way  around,  each  data  poim  not  on  a  slice  edge  is  connected  by  grid  lines  to  four 
of  its  neighbors  on  hs  own  dice,  plus 
the  point  above  or  below  h).  The 
lattice  divides  the  space  of  the  scalar 
fidd  into  cubes,  each  of  which  has  12 
grid  line  edges  and  8  data  point  comers. 

For  each  cube,  IC  determines  whether 
each  of  its  dght  comers  is  inside  or 
outside  the  level  surface.  This  information  b  used  to  infer  the  path  of  the  surfi^e.  In 
panicular,  if  some  comers  are  inside  and  some  outside  the  surfiice  fi^r  a  given  cube,  the 
surface  is  assumed  to  pass  through  the  cube.  Every  grid  line  shared  by  two  data  points  on 
opposite  rides  of  the  level  contributes  a  poim  on  the  surfiwe  where  it  cuts  the  grid  line. 

An  8-bit  binary  number  is  sufiidem  to  describe  every  possible  arrangemoit  of  inside  and 
outside  comers  and  therefore  every  posrible  inference  about  the  path  of  the  surfiue 
through  a  cube.  IC  uses  such  a  nunto  as  an  index  imo  a  table  of  triangular  surftce 
patches.  The  triangle  vertices  of  the  patches  all  lie  on  grid  lines  with  inside/outride  comer 
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pairs,  at  points  where  the  surftce  cuts  the  grid  lines.  The  exact  location  of  the  vertices  is 
determined  by  linear  interpolation  between  data  points  and  along  grid  lines. 

For  ead)  combination  of  scalar  field  and  levd,  1C  ends  up  with  a  list  of  triangles  that 
^iproximate  the  level  surfiice.  ICCanvas  then  uses  the  3D  graphics  functions  in  the  IC 
DLI.  to  present  interactive  views  of  this  surfiice  to  the  user.  Consistem  with  IC’s  design 
goals,  the  ICCarivas  (fi^lay  is  &st  and  ample.  Both  the  display  and  the  sur&ce  geometry 
can  be  saved  to  files  loaded  into  other  programs  that  provide  more  sophiaicated 
visualization  tools,  including  animation,  photorealistic  rendering  and  multiple 
semitransparent  sur&ces. 
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tfpBioai 

1C  File  Format 


This  q>pendix  discusses  the  binary  format  for  input  files  used  by  ICCanvas  and  the  IC 
DLL.  It  assumes  a  knowledge  of  programming  and  some  experience  with  reading  and 
writing  binary  files.  A  knowledge  of  the  C  language  and  of  3D  graphics  concepts  vtnll  be 
helpful. 

The  binary  format  of  simple  data  types  in  an  IC  data  file  is  based  on  the  internal  formats  of 
most  popular  C  compilers  running  on  Intel  80x86  PCs.  The  byte  order  places  the  most 
significant  byte  at  the  highest  linear  address,  an  order  sometimes  called  big-endian.  As 
u^  here,  the  C  types  short  and  Long  are  2-  and  4-byte  two’s  complement  imegers.  The 
C  type  float  is  a  4-byte  IEEE  floating-point  number.  The  header  structure  described 
below  is  packed,  meaning  that  no  extra  bytes  are  allowed  between  structure  members. 
The  file  format,  however,  aligns  all  data  in  tte  file  on  2-byte  boundaries. 

IC  data  files  con^  of  3D  density  data  preceded  by  a  118-byte  header.  The  data  values 
are  assumed  to  be  arranged  in  a  3D  lattice  or  structured  grid  (picture  the  steel  skeleton  of 
a  building).  Each  data  value  is  associated  with  a  point  in  space  where  the  grid  lines  cross. 
Data  poims  are  at  equal  intervals  for  a  ^ven  axis,  but  for  differem  axes  the  spacing  can  be 
different.  The  header,  defined  below,  gives  the  size  of  the  grid,  the  pacing  of  the  data 
points  along  each  axis,  and  other  information  about  the  data. 

typedef  struct  { 

unsigned  short  size; 

float  scale; 

float  offset; 

>  AXIST; 

typedef  struct  < 

char  comientC  40  3; 

char  reserved; 

unsigned  short  percentileC  9  3; 

long  creationjdate; 

float  simulati'on  tine; 

AXIST  axisC  3  3;” 

>  HEADER; 

comment  A  40'Character  NULL-terminated  string.  Intended  for  use  as  a 

name  or  description  field,  this  string  is  displayed  in  the  Level  and 
Progress  dialogs.  There  is  no  restriction  on  its  format. 

reserved  This  space  is  currently  ignored;  you  may  hide  what  you  want  here. 

percentile  An  array  containing  an  approximate  distribution  for  the  denrities. 

At  each  decile  interval  from  90%  to  10%,  the  array  gives  the  level 
at  which  that  percentage  of  the  non-zero  data  points  will  be 
enclosed  in  an  isosurfoce  at  that  level.  This  is  the  information 
displayed  by  the  Level  dialog.  Levds  are  floating-poiiit  values 
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written  in  the  same  16-bit  fonnat  described  below  for  the  density 
data. 

creation_date  This  is  a  creation  date  and  time  written  as  a  4-byte  integer  in  the 
manner  of  the  ANSI  function  timeO,  which  returns  the  number  of 
seconds  since  00;00;00  January  1,  1970.  It  is  meant  to  document 
when  the  original  data  was  created,  but  may  also  be  used  to  record 
when  the  IC  format  file  was  generated  from  the  data. 

siniulation_tiine  If  appropriate  for  the  data,  this  stores  the  time  in  floating-point 
seconds  fiom  the  time  origin.  For  simulations  that  produce  data  at 
a  sequence  of  time  points,  this  field  can  be  used  to  give  the  time 
coordinate.  The  simulation  time  is  displayed  in  the  Level  and 
Progress  dialogs. 

axis  An  array  of  three  structures  giving  the  dimensions  (size),  distance 

between  data  points  (scale),  and  offsets  from  the  origin.  The  scale 
and  ofibet  values  are  in  arbitrary  units.  Having  a  scale  for  each  axis 
allows  ICCanvas  to  anisotropically  scale  its  surfaces  for  cases 
where  the  volumetric  cells,  or  voxels,  from  which  the  data  have 
presumably  been  sampled  are  not  cubes.  The  offsets  allow 
ICCanvas  to  presove  spatial  information  in  cases  where  a 
simulation  has  tracked  a  moving  mass  to  keep  it  within  the 
bounding  box  of  its  space.  Both  scale  and  ofi^  are  taken  into 
accoum  when  a  surface  is  saved  with  the  Save  Surface  menu 
options. 

The  header  is  followed  by  the  actual  density  data.  Denrities  are  expressed  in  arbitrary 
units,  although  by  convention  they  have  been  in  grams  per  cubic  cemimeter.  Each  density 
value  is  a  floating-point  number  packed  into  16  bits  by  the  following  expression; 

X  <«  1e-16  ?  0  :  2048  *  loglOC  x  )  <*■  32768 

Packed  densities  may  range  fi'om  le-16  to  le+16,  with  a  precision  of  about  S  decimal 
digits  in  log  space.  The  special  value  0  is  reserved  for  “empty”  parts  of  the  volume.  The 
advantages  of  this  encoding  are  that  we  preserve  rdative  precision  over  a  large  range,  we 
store  only  the  precision  we  can  use,  and  we  are  able  to  treat  the  density  values  internally 
as  integers,  improving  performance  and  conserving  memory. 

Once  packed,  the  densities  are  run-length  encoded.  RLE  eliminates  redundancy  in  the 
data  array  by  replacing  runs  of  three  or  more  equal  values  with  a  count  followed  by  a 
single  instance  of  the  value.  Each  x-y  plane,  or  slice,  of  denrities  is  encoded  separately 
(runs  don’t  cross  slices)  and  is  preceded  in  the  file  by  a  two-byte  size  field  containing  the 
size  of  the  slice,  in  short  integers,  after  encoding. 

BASIC  and  C  source  code  for  producing  IC  format  data  files  is  included  on  the 
distribution  disk.  The  C  version  includes  several  machine-independent  functions  for 
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writing  the  binary  images  of  simple  daU  types,  as  well  as  a  function  called  PackQ  for 
performing  run*length  encoding.  A  simple  implementation  that  doesn’t  require  RLE  can 
be  written  on  the  basis  of  the  information  given  here.  When  writiiig  densities,  write  for 
each  slice  the  2-byte  size  1S1ZE*JSIZE41,  followed  by  a  single  2-byte  RLE  code 
ZS1ZE*JS1ZE-1,  followed  by  the  density  values  for  the  slice.  The  RLE  code  tells 
ICCanvas  to  treat  the  entire  slice  as  one  verbatim  block. 


Apromac 

ICMAKLBAS 


This  ^>peiidix  lists  a  BASIC  prognun  for  producing  data  files  in  IC  format.  To  run  the 
program,  type  qbaszc  at  the  DOS  prompt.  This  starts  the  BASIC  imetpreter  included 
with  versions  S.O  and  later  of  MS-DOS.  Open  ICMAKE.BAS,  the  text  of  the  program,  fi-om 
the  SOURCENicmakE  directory  on  the  distribution  disk,  then  press  F5  (Run).  A  sample 
input  file,  called  11EST.DAT,  is  included  on  the  distribution  disk. 


ICHAKE.BAS  . 

EW  22  Mar  94 

A  BASIC  prograa  for  creating  data  files  in  1C  foraat. 

This  prograa  deaonstrates  the  siaplest  approach  to  the  creation  of 
binary  1C  input  files  using  the  Microsoft  QBASIC  interpreter,  part 
of  MS-DOS  5.0. and  later.  ICMAKE.BAS  reads  a  text  file  containing 
30  density  data,  preceded  by  header  iteas  describing  the  data: 

naae  or  coaaent  text,  less  than  40  characters,  on  line  1 
a  siauLation  tiae  in  seconds,  or  0.0,  on  line  2 
for  each  axis:  an  array  diaension,  distance  between  data  points, 
and  origin  coordinate  on  lines  3,  4  and  5 
density  values  separated  by  coaaas 

It  then  writes  an  1C  binary  file.  Because  of  the  relative  slowness 
of  interpreted  BASIC,  run- length  encoding  is  not  done.  Each  slice 
is  instead  preceded  by  a  single  RLE  code  that  tells  1C  that  the 
entire  slice  is  in  the  file  verbatia. 

This  prograa  aay  be  tested  by  using  the  exaaple  text  file  called 
TEST.DAT  as  the  input  file. 


IC  header  definition 


TYPE  AxisT 
size 
scale 
offset 
END  TYPE 

TYPE  Header 
comment 
reserved 
percenti Le 
creatdate 
simtime 
xaxis 
yaxis 
zaxis 
END  TYPE 


AS 

INTEGER 

AS 

SINGLE 

AS 

SINGLE 

AS 

STRING  * 

40 

AS 

STRING  * 

22 

AS 

STRING  * 

18 

AS 

LONG 

AS 

SINGLE 

AS 

AxisT 

AS 

AxisT 

AS 

AxisT 

unsigned  shortC  9  3 
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*  .  procedure  decleratione 

DECLARE  SUB  GetHeader  (hdr  AS  Header) 

DECLARE  SUB  GetSLice  (nf  AS  INTEGER,  bufO  AS  INTEGER,  size  AS  INTEGER) 
DECLARE  SUB  PutHeader  (nf  AS  INTEGER,  hdr  AS  Header) 

DECLARE  SUB  PutSlice  (nf  AS  INTEGER,  bufO  AS  INTEGER,  size  AS  INTEGER) 
DECLARE  SUB  MakeDist  (hdr  AS  Header) 

DECLARE  SUB  SuaDist  (hdr  AS  Header,  bufO  AS  INTEGER,  sz  AS  INTEGER) 
DECLARE  SUB  ShowProgress  (i  AS  INTEGER,  iaax  AS  INTEGER) 

DECLARE  FUNCTION  LeveLl2D!  (i  AS  LONG) 

DECLARE  FUNCTION  LeveLD2IX  (d  AS  SINGLE) 


DIP  k  AS  INTEGER,  sz  AS  INTEGER 
DIM  hdr  AS  Header 
DIM  perc(1024)  AS  INTEGER 
DIM  nperc  AS  LONG 


'  .  program  begins  here 

GetHeader  hdr 

sz  •  hdr. xaxis. size  *  hdr.yaxis.si 
DIM  src(sz)  AS ‘INTEGER 
PutHeader  2,  hdr 

FOR  k  ■  1  TO  hdr. zaxis. size 
GetSLice  1,  sreO,  sz 
SunOist  hdr,  sreO,  sz 
PutSlice  2,  sreO,  sz 
ShowProgress  k,  hdr. zaxis. size 
NEXT  k 

MakeDist  hdr 
PutHeader  2,  hdr 

END 

'  .  of  ICMAKE 


*  slice  index  and  size 
'  IC  header 

'  density  distribution 
'  number  of  non -zero  data  points 


SUB  GetHeader  (hdr  AS  Header) 
DIM  s  AS  STRING 


PRINT  •' . •• 

PRINT  "  IC  Make" 

PRINT  " . " 

INPUT  "Input  File . ",  s 

OPEN  s  FOR  INPUT  AS  1 
INPUT  "Output  File . ",  s 


OPEN  s  FOR  BINARY  ACCESS  WRITE  AS  2 
INPUT  #1,  s 

hdr. comment  ■  LEFT$(s  +  CHRS(O),  40) 

INPUT  #1,  hdr.simtime 

INPUT  hdr. xaxis. size,  hdr. xaxis. scale,  hdr  .xaxis. offset 
INPUT  hdr. yaxis. size,  hdr. yaxis. scale,  hdr .yaxis. offset 

INPUT  #1,  hdr. zaxis. size,  hdr. zaxis. scale,  hdr. zaxis. offset 
END  SUB 


SUB  GetSLice  (nf  AS  INTEGER,  bufO  AS  INTEGER,  size  AS  INTEGER) 
DIM  i  AS  INTEGER 
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DIM  X  AS  SINGLE 

FOR  i  *1  TO  tilt 
INPUT  9r\i,  x 
buf(l)  >  LtvtL021(x) 
NEXT  i 
END  SUB 


FUNCTION  L«vcLD2IX  (d  AS  SINGLE) 

DIM  i  AS  LONG 

IF  d  <  1E-16  THEN 
1  -  0 
ELSE 

i  -  INT(2048«  *  (LOGCd)  /  2.30258S12497«>  *  32768#) 
END  IF 

LevelD2I  «  (i  AND  8HFFFF) 

END  FUNCTION 


SUB  MakcDist  (hdr  AS  Header) 

DIM  n  AS  LONG,  s  AS  LONG,  d  AS  LONG,  j  AS  LONG 
DIM  1  AS  INTEGER,  j2  AS  INTEGER 
SHARED  percO  AS  INTEGER 
SHARED  nperc  AS  LONG 

d  ■  nperc  \  10 
s  -  0 
i  -  8 

FOR  n  ■  2  TO  1024 
s  *  s  ♦  perc(n) 

DO  WHILE  s  >>  d 

j  ■  64  *  n  -  <64  *  <$  •  d))  /  perc(n) 
j2  >  j  AND  SHFFFF 

MlDSChdr. percentile,  i  *  2  1,  2)  >  MKIS(}2) 

s  *  s  -  d 
i  *  i  -  1 

IF  i  <  0  THEN  EXIT  DO 
LOOP 

IF  i  <  0  THEN  EXIT  FOR 
NEXT  n 
END  SUB 


SUB  PutHeeder  (nf  AS  INTEGER,  hdr  AS  Header) 
PUT  #nf,  1,  hdr 
END  SUB 


SUB  PutSlice  (nf  AS  INTEGER,  bufO  AS  INTEGER,  size  AS  INTEGER) 
DIM  rlecode  AS  INTEGER,  i  AS  INTEGER 

rlecode  *  size  •  1 
i  ■  size  +  1 
PUT  #nf,  ,  i 
PUT  #nf,  ,  rlecode 
FOR  i  *  1  TO  size 
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buf(i) 


PUT  «nf^  ^ 
NEXT  1 
END  SUB 


SUB  ShowProgress  (i  AS  INTEGER^  iMx  AS  INTEGER) 

CONST  F6GAUGE  «  178 
CONST  BKGAUGE  >  176 

IF  1  >  1  THEN  PRINT  STRINGS(21,  BKGAUGE); 

LOCATE  ^  1:  PRINT  STRINGS((21  *  i)  \  imat,  FGGAUGE); 
IF  (i  >  imax)  THEN  PRINT  :  PRINT  “Done." 

END  SUB 


SUB  SuaDist  (hdr  AS  Header,  bufO  AS  INTEGER,  sz  AS  INTEGER) 
DIM  i  AS  INTEGER,  j  AS  INTEGER 
SHARED  percO  AS  INTEGER 
SHARED  nperc  AS  LONG 

FOR  1  ■  1  TO  sz 

IF  buf(0  <>  0  THEN 

j  -  (buf(i)  AND  8HFFFF)  \  64 
perc(j)  «  perc(j)  ♦  1 
nperc  ■  nperc  +  1 
END  IF 
NEXT  i 
END  SUB 
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